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In [l[ Shao e£ a/, claim, based on low statistics simula- 
tions, that a model with majority rule coarsening exhibits 
in d = 2 a percolation transition in the universality class 
of invasion percolation with trapping (IPT). They report 
also that the system reaches its final state rapidly with 
no diverging time scale. Since the original configurations 
are random and thus in the ordinary percolation (OP) 
universality class, it seems unlikely that long range cor- 
relations could develop in a finite time that would change 
this. Indeed it was proved rigorously Q that similar 2-d 
models (called "dependent percolation" in belong to 
the OP universality class. 

Here we present high statistics (up to L = 2 14 , > 10 4 
realizations) on L x L square lattices and confirm that 
the phase transition is in the OP universality class, thus 
refuting a central tenet of [l[ . Initially each site i is ran- 
domly assigned one of two opinions (or spins): Ui = +1 
with probability /, otherwise <Xj = — 1. At each time 
step all sites arc updated in parallel. If at least three 
of their four neighbors disagree with them, they change 
their opinion, otherwise they keep it. As noted in [l| this 
leads quickly (within O(10) time steps) to a static state, 
except for sites that flip permanently with period 2. The 
critical probability f c where a cluster percolates de- 
pends slightly on how these flicker sites are treated (we 
treat them as "+", if <x; = +1 at even times), but the 
universal properties do not. 

We first determine f c by measuring the chance that a 
cluster in the final state percolates through lattices with 
open boundary conditions. Using finite size scaling Q, 
we obtain f c = 0.506425(20), in agreement with the less 
precise estimate of After that, we measure the dis- 
tribution of cluster sizes with a — +1 in final states ob- 
tained with helical b.c. for / « f c . Figure 1 confirms 
the above estimate of f c and shows that the data are 
excellently described by a power law P(s) ~ s~ T with 
the OP critical exponent t = 187/91 « 2.055, ruling 
out the IPT exponent r « 1.89. We see also deviations 
from this power law at small masses s, as small clus- 
ters are eliminated by the coarsening. This, together 
with using open boundary conditions and neglecting fi- 
nite size corrections, explains why r was underestimated 
in [l[. A data collapse of the r.h.s peaks in Fig.l gives 
Df = 1.895(15) as for OP, but in disagreement with IPT. 
Notice that the exponents obtained in [l[ strongly violate 
the hyperscaling relation r = d/Df + 1. 


Shao et al. claim that IPT is relevant because local 
clusters get trapped. The difference between OP and IPT 
is that clusters can grow both outwards and inwards (into 
empty holes) in OP, while they can only grow outward 
in IPT. In this respect the model of [l| is exactly as OP. 

We also simulated the process on random Erdos-Renyi 
networks. For small average degrees we confirm the claim 
of [l[ that the percolation transition is in the OP class. 
But for large average degrees we find an unexpected first 
order transition Q. 

As a model for opinion dynamics the model is of limited 
interest, since the dynamics leaves essentially unchanged 
all large clusters present in the initial state - except for 
clusters with hubs, in case of scale-free networks, who 
immediately adopt the majority opinion. 



FIG. 1. (Color online) Probability distribution of cluster sizes 
at / = 0.506425 for different L on a log-log scale. The straight 
lines represent power laws with exponent r = 2.055 (red, 
dashed) and r = 1.89 (green, dotted), corresponding to OP 
and IPT, respectively. Inset: Data for L = 4096 at different 
values of /, after multiplication with s 2 ' 055 . Small changes of 
/ give rise to deviations from power law behavior. At our es- 
timate of f c obtained independently from the spanning prob- 
ability, power law scaling with the OP exponent r = 2.055 
extending over three orders of magnitude. 

We thank a referee for pointing out the violation of 
hyperscaling in [l|. 
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